Quantum phase retrieval of a Rydberg wave packet using a half-cycle pulse 

J. Aim, D. N. Hutchinson, C. Rangan, and P. H. Bucksbaum 
Physics Department, University of Michigan, Ann Arbor, MI 48109-1120 

A terahertz half-cycle pulse was used to retrieve information stored as quantum phase in an 
A''-state Rydberg atom data register. The register was prepared as a wave packet with one state 
phase-reversed from the others (the "marked bit"). A half-cycle pulse then drove a significant 
portion of the electron probability into the flipped state via multimode interference. 

PACS numbers: 32.80.Rm, 03.67.-a, 42.30.Rx 
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An atomic Rydberg wave packet is an atom with one 
of its electrons in a coherent superposition of many high- 
lying orbitals. Wave packets can be viewed as data reg- 
isters for the information contained in the relative quan- 
tum phase of their constituent orbitals. One problem is 
that phase is difScult to detect. State selective detec- 
tion on an ensemble of identically prepared registers usu- 
ally only reveals the component amplitudes; the phases 
are hidden. Optical techniques have been employed to 
store and then efficiently retrieve phase information from 
atomic Rydberg wave packets [Q , following a suggestion 
that binary information stored as phase in quantum reg- 
isters could be searched more efficiently than the same 
classical information in classical data registers [EJ. This 
Letter reports a different technique for extracting phase 
information, based on the interaction between the wave 
packet and a broadband coherent electric field pulse. 

In the experiment reported in Ref. Q|, a shaped op- 
tical pulse carried binary information into the Rydberg 
atom, by creating an electron wave packet with one or 
more states phase- ffipped (binary 1) with respect to the 
others (binary 0). A second optical pulse created an- 
other wave packet (the decoder) which holographically 
interfered with the register wave packet [HQ] . This inter- 
ference converted the phase information to amplitudes 
by amplifying the flipped states and suppressing the un- 
flipped states. This method of information storage is 
limited because it depends on mutual coherence between 
a low-lying launch state and the Rydberg states. Fur- 
thermore, the decoding pulse requires that most of the 
electron's probability amplitude resides in the launch 
state, and only a small portion goes into Rydberg states. 
Therefore, the scheme demonstrated in the above exper- 
iment cannot be scaled to a very large data register g] . 

In the new data storage and retrieval experiment pre- 
sented in this paper, the launch state is no longer part 
of the data register, and all of the probability resides in 
the Rydberg states of the atom. The decoding process 
of amplifying the flipped bits is now performed by a tera- 
hertz domain half-cycle pulse (HCP) M. This broadband 
coherent far-infrared radiation directly couples together 
Rydberg levels. Information stored in the phase domain 
is converted to amplitude information by coherent redis- 
tribution induced by the HCP. 



HCP interactions with the Rydberg atoms have been 
studied previously, in HCP-ionization experiments and 
model calculations 0,§1, HCP-redistribution of energy 
eigenstates |9|] , and interactions with Rydberg wave pack- 
ets [|lO|,Q. In the present work, an HCP amplifies a 
phase-flipped orbital (the marked bit), suppressing the 
other unmarked states in the wave packet. 

The HCP locates the marked bit because its activity 
depends on the wave packet's location. The interac- 
tion is strongest when the wave packet is localized at the 
ion core, and this situation corresponds to a coherent 
sum of orbitals with equal phase Q . A Rydberg wave 
packet data register with one flipped bit is described by 
a superposition of such a localized wave packet, plus a 
phase-reversed orbital: 
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A relatively weak HCP interaction will destroy the lo- 
calized component easily, while, the spatially delocalized 
single marked orbital \nolm) is amplifled by a factor of 
approximately 2^ = 4. This is an example of the quan- 
tum diffusion and quantum ampliflcation of information, 
discussed by Grover Q . 

The database search experiment was performed in a 
beam of atomic Cs. The Rydberg data register was 
loaded using wave packet sculpting techniques that have 
been described previously [Q. Briefly, an optical pulse 
was sculpted with an acousto-optic Fourier filter |Q to 
excite a Rydberg data register with one marked bit in 
Cs pz-states with n=24-29. State-selective field ioniza- 
tion (SSFI) [n4| of the initial wave packet reveals only 
the amplitudes of the states, as shown in Fig. 1(a). 

An HCP was generated inside the atomic beam appara- 
tus by an ultrafast optical pulse illuminating a large aper- 
ture (1.5 X icm?) electrically biased GaAs wafer M. The 
far-infrared radiation was focused on the atomic beam 
by an off-axis parabolic mirror. Its peak field strength 
was calibrated by observing the threshold for ionization 
of Rydberg states m . The HCP delivered an impulse of 
approximately Qo = 0.0043?i/ao, which was enough to 
significantly mix s, p, and d states M. It was polarized 
parallel to the Rydberg state. 



The HCP interrogated the system at time r, convert- 
ing the wave packet to a simpler structure in which only 
the marked state remains. Figure 1 shows SSFI spectra 
for three different data register configurations, in which a 
different state in the register was phase-reversed as shown 
in Eq. (|l|) : the flipped bit was 25p, 26p, and 27p in panels 
(ii), (iii), and (iv), respectively. 

The quantum diffusion and marked state amplification 
can be calculated in the limit where the duration of the 
HCP is much less than the natural time scale for the 
phase evolution of the Rydberg wave packet, i.e. the Ke- 
pler orbit time, Ik — 27r/A_E — 27r(ri — Si)^, where n 
is the principal quantum number and Si is the quantum 
defect for the l^'^ angular momentum states. For atomic 
cesium at n = 26, tx ~ 2ps is five times the full width 
haff maximum (FWHM) of the HCP tfwhm ~ 440fs. In 
this impulsive limit, the HCP atomic interaction Hamil- 
tonian (in atomic units with e = rrie = h = I) 



H{t) = E{t) ■ r= -Q6{t) 
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where Q = j^_^ E{t')dt' is the total momentum trans- 
fer to the electron. The initial electron wave function 
\^i) acquires an extra phase along the direction of the 
impulse(z) 
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For an initial wave packet |^i) = ^„ an\nlm) , the final 
state wave packet is 
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where the transition matrix elements are given by 

/";™"' = {n'l'm'\e^-^\nlm). (5) 

For I = 1 and ni = 0, and Q directed along z, we find 

-il' - l)ju+iiQr)]\Rr,pir))6rn',o- (6) 

The phase of each transition f^i^ [Q) is determined 
by the parity change, Al = I' — I and the magnitude of 
Q. For example, p-p transitions are always real and p-d 
transitions are imaginary, while the sign changes with Q. 
The probability density in a final state \n'l'm') is given 

by 
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For atomic cesium levels around n w 26 excited by an 
impulse Q = Qo, the self-transition amplitude fl^^Q is n 



out of phase with transition amplitudes from other or- 
bitals /" Q . This creates constructive multi-mode inter- 
ference for the marked states, while transitions to other 
states experience destructive interference. The results 
of a calculation in the impulse limit are shown in Fig. 1, 
and show qualitative agreement. Deviations from the im- 
pulse limit were studied with a second calculation, where 
the Schrodingcr equation was integrated over a finite- 
duration model half-cycle pulse. We used an HCP model 
fitted to experimentally measured spectra [Q. The mo- 
mentum transferred by the HCP, Q, is the area under 
the curve. 
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The field-free basis states were calculated using a grid- 
based pseudopotential method |1^. The time-evolution 
was carried out by using a split-operator method in a re- 
stricted basis of essential states |lg] . This symmetrized- 
product propagator is unitary and correct through the 
second-order in the time step. We used a basis of 187 
states, that is, states with quantum numbers n = 21 — 31 
and i < 17{m ~ 0) and a time step of lOfs. The results 
are shown in Fig. 1. When a single initial state out of 
the wave packet is phase-flipped, we reproduce the result 
that this orbital survives the HCP with a significantly 
larger population than the rest, in good agreement with 
the experiment. 

The full dynamics of the data base are revealed by 
mapping the HCP redistribution as a function of the tar- 
get time. Since the Rydberg levels in the data register are 
non-degenerate, different configurations of flipped bits 
occur at times other than r. This two-dimensional phase 
space interference pattern, or "quantum carpet" p7| , is 
shown in Fig. 2, for data collected over several Kepler 
periods. The HCP consistently redistributed the wave 
packet into a single state when it consisted of only a sin- 
gle flipped orbital. The patterns of ridges in the quantum 
carpet show how the Kepler period changes with energy, 
and also display wave packet dispersion over 8 ps of the 
evolution. The lines are equivalent to classical ridges of 
a multi-mode interference pattern |17] ]. 

A quantitative measure of success in the database 
search problem is the reduction of the informational en- 
tropy of the system following the search algorithm [Q . 
Entropy S in this simple case of a single marked state is 
defined as 
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where Pi is the probability that the marked bit resides 
in orbital i. The information system consists of the A'^- 
state data register. When the data register is loaded, the 



information is encoded as the phase of each orbital. The 
orbital amplitudes carry no information. Since SSFI only 
measures amplitudes, the information is completely hid- 
den from view, and therefore the entropy of the system 
is maximum: S'initiai = -log(l/A^) = log TV. 

The interaction with the HCP converts the phase in- 
formation into amplitudes, which can be read out using 
SSFI. The rules of quantum measurement insure that 
only one state will be detected each time an atom is 
interrogated by SSFI. If the algorithm works perfectly, 
the marked orbital is revealed with unit efRciency, and 
the entropy drops to zero. Otherwise, unmarked orbitals 
may occasionally be detected by mistake. The entropy 
reduction can be calculated from the probability distri- 
bution, using Eq. (|9|). 

The states that are not included in the data register 
can rob some of the probability during the HCP oper- 
ation, and therefore contribute to the inefficiency of the 
algorithm. A lower bound on the amount of entropy con- 
tributed by leakage into the non-register states can be 
measured by including a reservoir state in the calcula- 
tion of the entropy. Table 1 shows the entropy reduction 
for each of the marked orbitals shown in Fig. 1. There is 
a significant entropy reduction, but the algorithm is far 
from perfect. Further progress could be made by shaping 
broadband terahertz radiation (which has been demon- 
strated using optical pulse shaping techniques [|l^) to 
focus the state redistribution. 

In conclusion, we have used a HCP to search for infor- 
mation stored as phase in a Rydberg wave packet. Phase 
information coded in a Rydberg quantum register was 
retrieved through the impulsive interaction of HCP. This 
interaction forced a constructive multi-mode interference 
only on a phase-ffipped state. Direct integration of the 
time-dependent Schrodinger equation as well as an im- 
pulsive model calculation showed good agreements with 
the experimental data. Beyond the problem of database 
searching, these results point toward the use of HCPs as 
general tools of analysis for unknown wave packets. This 
could be used in conjunction with optical techniques for 
wave packet tomography y. 
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TABLE I. Informational entropy calculations using Eq. H 
for experimental data Sexp, the impulse model calculation 
Simp, and the full calculation with the time-integration of the 
Schrodinger equation Sfuii. The initial entropy, S'initiai = 1.79 
for iV = 6. 
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FIG. 1. Experimental data and model calculations of 
the state-selective field ionization spectra of the ini- 
tial wave packet(i), and the final wave packets at 
r — 2.1ps(ii), 4.2ps(iii) and 4.7ps(iv). (a) Experimen- 
tal data (b) Impulse model calculation (c) Model calcu- 
lation includes the direct integration of the time depen- 
dent Schrodinger equation using the split-operator method. 
Both calculations assumed an initial Gaussian wave packet 
|Vt=o) = 0.5\2Ap) + \25p) + 1.2\2Qp) + \27p)+0.7\28p)+0.5\29p). 
The electron probability of the initial Gaussian wave packet 
was driven into the 25p(ii), 26p(iii), and 27p(iv) states re- 
spectively. The nearly degenerate p and d states and the s 
and hydrogen-like states are not resolved in the experimental 

data. 
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FIG. 2. Experiment(a) and theory(b) of the multi-mode 
interference pattern of Rydberg wave packets forced by the 
half-cycle pulse interaction. Guide lines follow the phase gra- 
dient maxima that appear around the half multiples of the 
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(fc -I- l/2)fK(n) with fc = 0, 1 and 2. 
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